Detailed Topography of the Fermi Surface of Sr 2 Ru04 
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We apply a novel analysis of the field and angle dependence of the quantum-oscillatory amplitudes in 
the unconventional superconductor Sr2Ru04 to map its Fermi surface (FS) in unprecedented detail, 
and to obtain previously inaccessible information on the band dispersion. The three quasi-2D FS 
sheets not only exhibit very diverse magnitudes of warping, but also entirely different dominant 
warping symmetries. We use the data to reassess recent results on c-axis transport phenomena. 

PACS numbers: 71.18.+y, 71.27.+a, 74.25. Jb 
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The layered perovskite oxide Sr2Ru04 has attracted 
considerable experimental and theoretical attention since 
the discovery of superconductivity in this compound five 
years ago M. Fermi liquid behaviour of several bulk 
transport and thermodynamic properties was observed 
in early work and the existence of mass enhanced 

fermionic quasiparticles was demonstrated explicitly by 
the observation of quantum oscillations in the magnetiza- 
tion (de Haas- van Alphcn or dHvA effect) and resistivity 
Quantitative similarities between the Fermi liquid in 
He and that in Si~2Ru04 hint at the possibility of p-wave 
superconducting pairing M. Evidence supporting such a 
scenario has come from the existence of a very strong im- 
purity effect |p| , a temperature independent Knight shift 
into the superconducting state a muon spin rotation 
study indicating broken time reversal symmetry j7j and a 
number of other experiments |^||. Taken together, these 
favour spin triplet superconductivity with a p-wave vec- 
tor order parameter and a nodeless energy gap. 

Sr2RuC>4 appears to be an ideal material in which to in- 
vestigate unconventional superconductivity in real depth: 
of all known compounds exhibiting this phenomenon, 
Si'2Ru04 offers the best prospects of a complete under- 
standing of the normal state properties within standard 
Fermi liquid theory |l0| . This would provide a solid foun- 
dation for all theoretical models. It has become clear, 
however, that further progress will require very detailed 
knowledge of the electronic structure of S^RuC^. For ex- 
ample, one of the most successful current theories relies 
on the assumption that the FS consists of two sheets de- 
rived from bands with strong Ru d xz ,yz character and one 
with strong Ru d xy character Q. These are supposed 
to form weakly coupled subsystems with very different 
pairing interactions This theory has successfully 

predicted non-hexagonal vortex lattice structures Jl3|,[l4| , 
but it is less clear whether, in its simplest form, it will pro- 
vide a satisfactory explanation for recent measurements 
of the temperature dependence of the density of normal 
excitations in the superconducting state [fl5| . 

Of central importance to the understanding of 



Sr2RuC>4 is the origin of the quasiparticle mass enhance- 
ment and how it relates to magnetic fluctuations. Re- 
cent observations of cyclotron resonances jl(| give the 
promise of separating the various contributions to the 
enhancement, but identifying the type of resonance and 
the extent to which electron interactions are affecting the 
observed masses requires more detailed knowledge of the 
Fermi surface than has been available to date. Clues to 
the magnetic fluctuation spectrum have come from nu- 
clear magnetic resonance p7| and neutron scattering ex- 
periments Q] and from calculations [[l9| . Both ferro- and 
antifcrromagnetic fluctuations appear to be present, the 
latter due to nesting of the FS. Angular dependent mag- 
netoresistance oscillations (AMRO) can give information 
about the in-plane topography of the FS and the extent 
to which it is nested. However, with three bands cross- 
ing the Fermi level, the interpretation of AMRO data on 
Sr2Ru04 pp| | has been somewhat ambiguous. 

Progress on all the issues discussed above requires 
high resolution, sheet-by-sheet knowledge of the FS of 
Sr2Ru04. As shown in previous studies ||, the dHvA 
effect is ideally suited to this, as data from individual FS 
sheets can be identified without ambiguity. For this rea- 
son, we have performed a comprehensive angular dHvA 
study in Sr2Ru04. Full analysis of the data required ex- 
tension and generalisation of previously reported theoret- 
ical treatments [^TJ of dHvA in nearly 2D materials. As 
a result, we present an unprecedentedly detailed picture 
of the warping of each FS sheet. A series of recent mea- 
surements of c-axis transport phenomena are discussed in 
light of the full experimentally determined dispersion. 

Quantum-oscillatory effects in a crystal arise from the 
quantization of the cyclotron motion of the charge carri- 
ers in a magnetic field B. For three-dimensional metals, 
only the extremal cyclotron orbits in fc-space lead to a 
macroscopic magnetization, and the quantitative treat- 
ment has been known for decades |^2| . 

For a quasi-2D metal, the FS consists of weakly corru- 
gated cylinders. While such weak distortions have little 
effect on the cross-sectional areas which determine the 
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FIG. 1. Example of a dHvA field sweep on Sr 2 Ru0 4 , for 6 = 9.8° off the c-axis in the [001] -> [110] rotation study. The 
vertical axis is the pick-up signal (in arbitrary units) at the second harmonic of the excitation frequency. At high fields, beats 
in the /3-oscillations are visible, as indicated by arrows. 



dHvA frequency, they still affect the interference of the 
magnetization contributions of different parts of the FS 
and therefore lead to a characteristic amplitude reduc- 
tion. Conversely, as we will show, analysis of the experi- 
mental dHvA amplitude behaviour can reveal fine details 
of the topography of the underlying Fermi cylinders. 

In the most basic case, a simple corrugation of the 
Fermi cylinder leads to a beating pattern in the magneti- 
zation. Analysis of the beats for on-axis fields gives some 
information about the magnitude of the warping, but fur- 
ther conclusions have to rely on assumptions about the 
precise form of the corrugation ||] . At the next level of ap- 
proximation, the FS dispersion can be determined within 
the traditional scope of the Yamaji effect [ pT|j2^] . A pre- 
liminary attempt to extract information on Sr 2 Ru04 in 
this way, however, has not achieved agreement between 
the data and the predictions of the simple model [p4f . 
We show that a much more extensive treatment is needed 
that (a) considers Fermi cylinder corrugation of arbitrary 
shape and (b) if necessary, goes beyond the extremal orbit 
approximation. Also, to extract meaningful information 
from experiments, one needs data of much higher quality 
than has been available to date. 

In the following, we will briefly present the results of a 
full quantitative treatment of the oscillatory magnetiza- 
tion for a Fermi cylinder that is warped arbitrarily but 
still compliant with the Brillouin zone (BZ) symmetry of 
Sr2RuC>4; details will be presented elsewhere p5[ . It is 
convenient to parameterize the corrugation of the cylin- 
der through an expansion of the local Fermi wavevector, 
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(see Table | for illustration). Here, k = ck z /2 where c is 
the height of the body-centered tetragonal unit cell, and 
<j) is the azimuthal angle of k in the (k x , fcj,)-plane. The (3- 
and 7-cylinders are centered in the BZ; symmetry allows 
nonzero fc M „ only for \i divisible by 4. The a-cylinder runs 
along the corners of the BZ and has nonzero fc M „ only for 
v even and /x divisible by 4, or for v odd and \x mod 4 = 2. 
The average Fermi wavevector is given by fcooj which is 
assumed to be much larger than the higher-order k^ v . 
One also has to consider the effect of the magnetic field: 



spin-splitting drives the spin-up and spin-down surfaces 
apart; the parameters k^ can be taken to expand weakly 
and linearly in the field as in kj} v — k^ u ± x^B. Ordi- 
nary spin-splitting is described by xoo alone, while higher- 
order contributions would correspond to the underlying 
electronic band structure being flatter at some points on 
the Fermi surface than at others. Indeed, this anoma- 
lous spin-splitting will prove essential for describing the 
a-sheet in Sr2Ru04. 

If a magnetic field is applied at polar and azimuthal 
angles 8 and </>, the Fermi surface cross-sectional area 
perpendicular to the field which cuts the cylinder axis at 
k is given by the Bessel function expansion 
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which is a generalization of Yamaji's earlier treatment 
plfl ; here, kf = cfcoo/2. The total quantum oscillatory 
magnetization now arises from the interference of the in- 
dividual contributions of the cylinder cross-sections; at 
constant chemical potential, the fundamental component 
of the oscillations is given by 
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Eqs. (||) and (||) describe oscillations at the undistorted 
frequency Hk 2 /2ecos6, with a characteristic amplitude 
modulation induced by the interference. One can thus 
infer the warping parameters k^ v by modeling experi- 
mentally obtained amplitude data (26) . 

We have performed a thorough dHvA rotation study 
in the [001] — > [110] plane on a high-quality crystal of 
Sr2Ru04 with T c > 1.3 K. The experiments were car- 
ried out on a low noise superconducting magnet sys- 
tem in field sweeps from 16 T to 5 T, at temperatures 
of 50 mK. A modulation field of 5.4 mT amplitude was 
applied to the grounded sample, and the second harmonic 
of the voltage induced at a pick-up coil around the sam- 
ple was recorded, essentially measuring d 2 M/dB 2 , with 
well-established extra contributions from the field modu- 
lation, impurity scattering, and thermal smearing B2|. 
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FIG. 3. Angular dependence of the beat frequency of the 
/3-oscillations in Sr 2 Ru0 4 (cf. Fig. |) in the [001] -> [110] rota- 
tion study (full circles). AF(0°) fixes k 01 at 3.9 x 10 7 m" 1 ; the 
data are compared against the predicted behaviour for cylin- 
drically symmetric warping (dashed line), for zero dispersion 
along the BZ axes (dotted line), and for fco4 as tabulated in 
Table § (solid line). 

A typical signal trace, demonstrating the high quality 
of our data, can be viewed in Fig. 0. Overall, 35 such 
sweeps were performed, at angular intervals of 2°. For 
each of the dHvA frequencies corresponding to the three 
FS sheets, we have extracted the local dHvA amplitude 
through filtering in the Fourier domain. The dHvA am- 
plitude can then be visualized versus magnitude and di- 
rection of the field, as shown in the density plot in Fig. || 
for the a-sheet. We have also performed similar data 
analysis for a second extensive rotation study of short 
high-field (18 T to 15 T) sweeps in the [001] -> [100] plane; 
these runs used a different sample. We now turn to the 
results of the analysis for all three FS sheets, where Ta- 
ble | presents the deduced values for the . 

a-Sheet — The most striking features of the data in 
Fig. are qualitative differences with similar data in the 
[001] -> [100] plane (Ref. @ and the present study), and 
the absence of spin-zeroes that should be visible as ver- 
tical black lines. We are able to account for both effects, 
and produce near-perfect agreement with experiments as 
seen by comparing the two panels of Fig. |^, using the k^ v 
as given in Table [|. The dominant fc2i-term affects the 
dHvA amplitude for <fr = 45° but has no effect for fields in 



FIG. 2. Density plot of the dHvA 
amplitude (in arbitrary units) of 
the a- frequency in Sr2Ru04, in 
the experimental [001] — > [110] ro- 
tation study (left), and compar- 
ison with the theoretical simula- 
tion, using the warping parame- 
ters in Table | (right). The the- 
oretical calculation incorporates 
experimental effects such as the 
field modulation amplitude char- 
acteristic. 

the [001] — > [100] plane. The absence of spin-zeroes arises 
from the presence of a finite X21 — ~5 x 10 5 m~ 1 T _1 in 
addition to xoo = 10-4 x 10 5 m _1 T _1 . The same param- 
eters account equally well for data from [001] — > [100] 
rotations and for the amplitude of the second harmonic 
for both rotation directions. It should be noted that the 
warping of the a-cylinder is so weak — at some angles 
and fields it is smaller than the Landau level spacing — 
that the success of the model requires the use of our exact 
treatment beyond the extremal orbit approximation. 

13-Sheet — The warping is comparatively large, and 
the extremal orbit approximation is valid over most of 
the angular range. The Yamaji angle of 32° and the vari- 
ation of the (relatively fast) beating frequency AF with 
9 (Fig. ||) reveal that the dominant warping parameters 
are ferji and ^41, as tabulated in Table ||. They have op- 
posite signs, so the c-axis dispersion is largest along the 
zone diagonals. It is difficult to extract meaningful in- 
formation on the higher-order terms, but we believe that 
the data set an upper bound for a double warping contri- 
bution of |fc 02 | < 10 7 m- 1 . The spin-splitting behaviour 
is intricate, and while it is certain that higher-order Xp,v 
are needed to explain the data, it is impossible at this 
stage to extract them without ambiguity. 

j-Sheet — Due to stronger impurity damping, the 7- 
signal is only observable at fields of more than 13 T. 
Along [001] -> [110], its amplitude peaks at 9 = ±13.7°, 
implying that the dominant corrugation is double warp- 
ing, i.e. fcrj2 3> fcoi- We obtained a rough estimate of 
its strength from the sharpness of this amplitude maxi- 
mum — it is difficult to assess fcce from an on-axis beat- 
ing pattern, as that cannot be established in the short 
field range over which 7-oscillations are visible. For the 
[001] — ► [100] rotation, the amplitude maximum occurs 
at 6 = 14.6°. The deviation of the two measured (3-values 
from each other and from the simple Yamaji predicition 
of 14.1° yields whose sign implies that the c-axis dis- 
persion is largest along the zone axes. At present, it is 
not possible to extract reliable information on the spin- 
splitting parameters. 

To the order of expansion given here, it is possible to 
calculate the contribution of each FS to the c-axis con- 
ductivity. The /3-sheet dominates with a 86% share, com- 
pared to 8% for the a- and 6% for the 7-sheet. This 
provides new insight into recent AMRO experiments by 
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TABLE I. Summary of the extracted warping parameters 
k uv (in units of 10 7 m _1 ) of the three FS sheets of S^RuCU 
[B6|. Entries symbolized by a long dash are forbidden by the 
BZ symmetry. The first row visualizes the warping for the 
different values of and v, on top of a large fcoo- 
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Ohmichi et al. |20| as it strongly suggests that the AMRO 
signal originates predominantly from the /3-sheet (rather 
than the a-sheet as had previously been assumed). The 
AMRO data then fix the "squareness" parameter of the 
/3-cylinder as k± — 5.3 x 10 8 m _1 , giving quantitative 
information about the FS nesting on that sheet. 

Our work also helps to clarify the interpretation of ab- 
sorption spectra in cyclotron resonance |l6| , which have 
recently been used to assess quasiparticle mass enhance- 
ments in Sr2Ru04. The "periodic orbit resonance" ge- 
ometry used in those experiments assesses modulations in 
the c-axis Fermi velocity. The strongest signals from such 
resonances in their simplest form should again stem from 
the /3-sheet, and the unusual warping symmetry would 
lead to dominant resonances at 4w c , 8lo c , . . . for the (3- 
and 7-orbits, and at 2w c , 4w c , . . . for the a-orbit. 

High-precision dHvA itself can provide normally inac- 
cesible information on spin-dependent many-body effects, 
by measuring both the specific heat and the spin suscep- 
tibility mass enhancements. For the a-shect, we have 
m*/m = 3.4 and m* usc /m = 4.1, the latter obtained 
from spin-splitting analysis assuming g ~ 2. We would 
expect the ratio m* usc /m* to diverge at a ferromagnetic 
quantum critical point; the small ratio here suggests that, 
at least for the a-sheet, the paramagnetic susceptibility 
enhancement is matched by specific heat contributions 
from phonons or large-g spin fluctuations [|l0| , p~8|| . 

Finally, an intriguing feature of the present study is 
the qualitative difference between the dominant warping 
symmetries observed for the three FS sheets of Sr2Ru04. 
Detailed comparison with the results of band structure 
calculations would be very informative to test the ac- 
curacy of these computations for the weak out-of-plane 
dispersions in quasi-2D metals in general. For Si~2Ru04 
in particular, this would also be a check on the commonly 
made assignment of the Ru d XZjyz orbital character to the 
a- and /3-sheets, and the d xy character to the 7-sheet, 
which lies at the foundation of theories of orbital depen- 
dent superconductivity. 

In summary, the full analysis of angle-dependent dHvA 
amplitude data has emerged as an extremely powerful 
tool to determine the exact topography of the FS in quasi- 
2D metals. We have been able to extract quantitative in- 



formation on the corrugation of all three Fermi cylinders 
of Sr2Ru04. The single warping of the /3-sheet provides 
most of the c-axis dispersion, while the ellipsoidal warp- 
ing of the a- and the double warping of the 7-sheet are 
less significant. 
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